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ABSTRACT 

Aims. We investigate the effects of magnetic fields on neutrino production in microquasars. 

Methods. We calculate the steady particle distributions for the pions and muons generated in pj and pp interactions in 
the jet taking the effects of all energy losses into account. 

Results. The obtained neutrino emission is significantly modified due to the synchrotron losses suffered by secondary 
pions and muons. 

Conclusions. The estimates made for neutrino fluxes arriving on the Earth imply that detection of high-energy neutrinos 
from the vicinity of the compact object can be difficult. However, in the case of windy microquasars, the interaction 
of energetic protons in the jet with matter of dense clumps of the wind could produce detectable neutrinos. This is 
because the pions and muons at larger distances from the compact object will not be affected by synchrotron losses. 

Key words. X-rays: binaries - neutrinos - radiation mechanisms: non-thermal 



1. Introduction 

Microquasars, the X-ray binary systems with non-thermal 
jets, are considered important candidate sources of high- 
energy neutrinos (Waxman & Levinson 2001). The recent 
detection of TeV gamma rays reveals that these objects 
are capable of accelerating particles to very high ener- 
' gies (Aharonian et al. 2005, Albert et al. 2006, Albert et 
al. 2007). The models that predict both gamma ray and 
neutrino emission are based on interactions of relativis- 
tic protons in the jet with cold protons of a dense wind 
from a high-mass stellar companion (Romero et al. 2003, 
Christiansen et al. 2006), with secondary synchrotron emis- 
sion in the jet itself (Romero & Vila 2008), and with cold 
protons in a heavy jet (Reynoso et al. 2008). A usual as- 
sumption made in these models is equipartition between the 
magnetic energy and the kinetic energy in the jets, which 
leads to large magnetic fields. In this work, we analyze the 
effects caused by the presence of such strong magnetic fields 
on the spectra of secondary particles that decay to neutri- 
nos. 

The outline of this work is as follows. In the next sec- 
tion we briefly discuss the basics of hadronic models for 
microquasars, and in Sect. [3] we deal with the acceleration 
and cooling mechanisms relevant to the primary relativistic 
particles in the jet. In Sect. 01 we analyze the effects of the 
magnetic field on the spectra of secondary pions, muons, 
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and neutrinos. In Sect. [5j we discuss the neutrino produc- 
tion through interactions between the jet and clumps of 
the stellar wind in high-mass microquasars. The last two 
sections include a discussion of the results and a summary. 

2. Basics of hadronic models of microquasars 

In these models, an accretion disk is present around the 
compact object, and a fraction of the accreted material is 
expelled in two oppositely directed jets (Falcke & Biermann 
1995). We assume conical jets with a half-opening angle £ 
and radius r{z) = zotan£, where the injection point is at 
a distance z$ from the compact object. A sketch of a high- 
mass microquasar is shown in Fig.[TJ where the star presents 
a wind with a clumply structure. In the case of a low-mass 
microquasar, there is no significant stellar wind. 

The kinetic luminosity of the jet, L^, implies a kinetic 
energy density of 



where Vh is the bulk velocity of the jet particles. Following 
the jet-accretion coupling hypothesis, we assume that 
around 10% of the Eddington luminosity goes into the jet 
(Kording et al. 2006). We adopt L k = 10 38 erg s" 1 for a 
10 Mq black hole. Equipartition then implies a magnetic 
energy density /O m ag = pk, and hence a magnetic field (e.g. 
Bosch-Ramon et al. 2006) 

B{z) = y/8irp k (z). (2) 
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Fig. 1. Schematic view of a high-mass microquasar. 



We consider that a fraction of the kinetic power in the 
jet is carried by relativistic primary protons and electrons, 
L lc i — L p + L e . The relation between the proton and elec- 
tron power is given by a certain parameter a in such a 
way that L p = a L e . This parameter is unknown, although 
there are reasons to think that a > 1. We consider the 
cases of a = 1 for equal proton and electron luminosities 
and a — 100 for a proton-dominated jet. 



3. Hadronic processes at the base of the jet 

In the one-zone approximation (Khangulyan et al. 2007), we 
assume that shock acceleration of the particles takes place 
in the jet at distances from zq to z max = 5zrj from the com- 
pact object. The injection rate is assumed to be a power law 
in the particle energy N'(E') = K E'~ 2 (GeY~ 1 cm- 3 ). The 
corresponding current can be written as J'(E') — cN'(E') 
in the reference frame co-moving with the jet particles. The 
conservation of the number of particles is satisfied if the 
current evolves with z as (see Ghisellini et al. 1985) 



J'(E', z) = K c (^) 2 S'-^GeV-VW- 2 ). 



(3) 



The continuity equation in the case of no time depen- 
dence and in the absence of sinks, implies that the injection 
or source function of particles must satisfy 

Q'(E',z) = V • J'(E',z) z, 

so that we have 

Q(E',z) = g (^) 3 ^'- 2 (GeV- 1 cm- 3 s- 1 ). (4) 

In the observer reference frame, whose line of sight 
makes an angle with the jet, since E' — T\,(E — 
(3 h cos 9y/E 2 -m 2 c 4 ), dV'/dV = T h , and 



dE' _ /3 h E cos 6 

~d~E ~ b ~ V^ 2 -™^' 



(5) 





Parameter 


Value 




Lk: jet power 


10 38 erg s" 1 


t 


q TC \: jet's content of relativistic particles 


0.1 




a: hadron-to-lepton ratio 


1, 100 


• 


zo : jet's launching point 


10 8 cm 




z mal : extent of acceleration region 


5z 


• 


Tb: jet's bulk Lorentz factor 


1.25 




£ : jet's half-opening angle 


1.5, 5° 


• 


6: viewing angle 


30° 




r): acceleration efficiency 


0.1 



it follows that 



Q(E,z) - Q (^7) 3 IV 1 (tf-zW^-mVcos 

(3b E cos 9 



r b - 



V^ 2 



where Tb is the bulk Lorentz factor of the jet. The normal- 
ization constant Qq is obtained by specifying the power in 
relativistic particles: 



dE e pE e pQ e p(E e p, z). 



The minimum energies are Eg 111 ' — 1 MeV and E i 
1.2 GeV, and the maximum energies will be obtained in the 
next section by equating the acceleration rate to the energy 
loss rate. The parameters of our model are summarized in 
Tabled! 



(7) 



n(min) 
P 



3.1. Accelerating and cooling rates: maximum energies 

The rate of acceleration of the particles to an energy E, 
t^ c = E~ 1 dE/dt, is given by 



t: 



ri- 



ce B 



En 



(8) 



where we consider r) = 0.1 for the acceleration efficiency. 
This corresponds to the case of an efficient accelerator, as 
expected at the base of the jet where shocks are mildly 
relativistic; see, e.g., Begelman et al. (1990). 

Charged particles of mass m and energy E — 7 mc 2 will 
emit synchrotron radiation at a rate 



r 1 = 

sync 



4 /m e y 
3 V m / 



cr T B 2 
m e c 8ir 



7- 



(9) 



In the jet at a distance z from the compact object, the 
density of cold particles is 



/ s (l-?rel) T 

1 rripC^irr^Vb 



(10) 



The rate of pp collisions of the relativistic protons with the 
cold ones is then given by 



t£ = n(z)c*<W(E p )K„, 



(11) 
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where the inelasticity coefficient is K pp w 1/2 and the cor- 
responding cross section for inelastic pp interactions can be 
approximated by (Kelner et al. 2006) 



°pp el) ( E P ) = (34.3 + 1.88L + 0.25L 2 ) 



Eth 



x 10 



- 2 w, 



(12) 



where L = ln(-E p /1000 GeV) and E th = 1.2 GeV. Because 
the jet is expanding with a lateral velocity (t> D tan£) the 
adiabatic cooling rate is (Bosch-Ramon et al. 2006) 



ad 3 z ' 



(13) 



We estimate the maximum energies achieved by the parti- 
cles by equating ^(£( roax )) = ^(E^^). In the case 



of electrons, we assume t 



loss 



t S y n + t ad , and for protons 



ti„ M = t„}„ + + t p £. At the base of the jet, for a = 100, 



loss "syn 

we obtain 



pp 



7 GeV 



(14) 



£> ax) (z ) 
and 

(z ) ps 10 7 GeV. (15) 

We show in Fig.[3]the above rates for electrons and protons, 
as well as the cooling rates due to pj and IC interactions 
that arise due to the photons from the synchrotron emission 
(see next section). 

3.2. Proton and electron distributions 

In the one-zone approximation the particle distribution in- 
dependent of time, i.e., in a steady state, can be obtained 
as the solution of the following transport equation: 



dN(E,z)b(E,z) 
dE 

where b(E, z) = 



+ t- s l(z) N(E,z)=Q(E,z), 
-Et^jE,z) and 



(16) 



(17) 



is the escape rate. 

The corresponding solution is 



N(E, z) 



1 



HE)\ 



dE'Q(E', z) 
xexp{-t£(z)T(E i E')}, (18) 



with 



r(E,E') 



E' 



dE" 



We notice that here the effect of particle acceleration is 
included through the injection function which depends on 
the energy with a power law in the frame co-moving with 
the bulk of the jet. 

We show the obtained distributions as a function of en- 
ergy and z in Figs. [3] and [4] for protons and electrons, re- 
spectively. It can be seen from the latter figure that if a = 1 
the number of electrons is higher than for a = 100, as ex- 
pected. For protons, we show the case of a = 100, but it 
does not differ significantly if a = 1. 




Fig. 3. Proton distribution as a function of energy and dis- 
tance to the compact object. 




8.2 \CM 



Fig. 4. Electron distribution as a function of energy and 
distance to the compact object. The cases of a = 1 and 
a = 100 are shown in gray and transparent surfaces, re- 
spectively. 

3.3. Synchrotron radiation 

Both the protons and electrons will emit synchrotron radi- 
ation. The power radiated by a single particle of energy E 
and pitch angle a is (e.g. Blumenthal & Gould 1970) 



P syn (E 7 ,E, z,a) 



V3e 3 B(z) Ey 
ATrmc 2 h E CI J E ^/ E , 



d(K 5/3 (t), (19) 



where K 5 / 3 (£) is the modified Bessel function of order 5/3 
and 

3heB(z)sma ( E 
mc 2 



E c 



The power per unit energy of the synchrotron photons is 



r( e >P) 



P syn N e}P {E, z), 



(20) 
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Fig. 2. Accelerating and cooling rates for protons (left panels) and for electrons (right panels) at the base of the jets. The 
top panels correspond to a — 100 and the bottom ones to a = 1. Plots shown: (solid lines), adiabatic cooling rates (dotted 
lines), synchrotron cooling rates (long-dashed lines), pp cooling rates (short-dashed lines), p-f cooling rates (dash-dotted 
lines, left panels), and IC cooling rates (dash-dotted lines, right panels). 



and the total luminosity can be obtained by integrating in 
the volume of the region of acceleration 



4y'i ) (E J )= [ d 3 rE,e^. (21) 
Jv 



The results for synchrotron radiation of protons and elec- 
trons are shown in Fig. [5] with a = 1 in the right panel and 
a = 100 in the left panel. 



3.4. Inverse Compton and pj interactions 

The synchrotron photons will, in turn, serve as targets for 
electrons and protons themselves. Locally, the correspond- 
ing radiation density can be expressed as 

N ph (e,z)^ £ -^ r -^-(GeV- 1 cm- 3 ). (22) 
e c 

Electrons will interact by Inverse-Compton scatterings 
at a rate 

t^{E,z) = \^ le , (23) 
3 m e c 
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Fig. 5. Synchrotron luminosity emitted by protons (solid lines) and by electrons (dashed lines). Black lines correspond 
to £ = 5° and green lines to £ = 1.5°. The cases of a = 1 and a — 100 are shown in the right and left panels, respectively. 



where 



e N ph (e)de 



But before decaying, pions and muons may interact, losing 
energy according to the processes discussed in the previous 
section. In the case of pions, 



is the corresponding energy density in photons. As protons 
interact with synchrotron photons, they lose energy due to 
photopion production at a rate 



, dE 
b^{E,z) = — 



ph 

2e7 p 



x / de'o${e')K${e')e'. (24) 



E(t S yn + ^ad + ^np ^7)- 

For the irp interactions we consider 
1, , _ n(z) ca^p Cl) {E p ) 

with a 



(29) 



(30) 



rr;™'(E) based on the proton being formed 



Here, e t h = 150 MeV and we adopt the cross section 
(Atoyan & Dermer 2003, see also Kelner & Aharonian 2008) 



'trp\^) ~ 3°pp 

by three valence quarks, while the pion is formed by two 
(Gaisser 1990). As for the 7r7 interactions, we estimate a 
cooling rate using expression (|24|) with the replacement 

For muons, we have 



Up-y 



(2/3)0$. 



4; ] = 0(c' - 200 MeV) 6(500 MeV- e') 3.4 x 10- 28 cm 2 z) = ~E(t^ n + tj + t lc L ). (31) 

In Fig.[S]we show the different rates corresponding to z = zq 
for pions in the left panels and for muons in the right panels. 
The cases with a = 100 are shown in the upper panels and 
the cases with a = 1 in the lower panels. We have included 
the rate of decay and escape as 

t-^(E,z)=t-l(z)+t^ c (E), (32) 

2.6 x 10~ 8 7 Tr ]^ 1 (s _1 ) for pions and t^ c = 
7 Al ]~ 1 (s~ 1 ) for muons. 



+ 9(e' - 500 MeV) 1.2 x 10" 28 cm 2 , (25) 
and the inelasticity as 

K$ = 6(e' - 200 MeV) 6(500 MeV - e') 0.2 

+ 6(e' - 500 MeV) 0.6. (26) 



The obtained IC and P7 cooling rates are shown in Fig. 
[2] It can be seen from this plot that the dominant mecha- 
nisms for energy loss are those discussed at the beginning 
of this section. 



where *dec 



4. Magnetic effects on secondary particles 

The primary relativistic protons will produce pions through 
inelastic interactions with matter and radiation. Pions will 
decay to muons and neutrinos, and muons will also decay, 
giving neutrinos and electrons: 

tt~ — > — > eTv^VeD^ (27) 
7T + — * /i+z^ — * e + v tl v e v tJ _. (28) 



4.1. Pion injection 

The injection function of pions produced by pp interactions 
is given by 



Q<r\E,z) = n{z)c 



dX ST f E 

—Np -,z 

X \ X 



(max) 
P 

xi^'UJJ <#->(§) (33) 
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Fig. 6. Cooling rates for pions (left panels) and for muons (right panels) at the base of the jets. The top panels correspond 
to a = 100 and the bottom ones to a = 1. Plots shown: adiabatic cooling rates (dotted lines), synchrotron cooling rates 
(long-dashed lines), wp cooling rates (short-dashed lines), irj cooling rates (dash-dotted lines, left panels), and IC cooling 
rates (dash-dotted lines, right panels). The decay plus escape rates are also shown (solid lines) 



where The injection function for charged pions from P7 inter- 

actions is 

E\ „ „ , / 1 -x a 



F (pp) (V Ej = 4aB 7T x a ~ 1 ( 



1 + rx a (l - x a ) / 



l-x a 1 + rx a (l - x a ) J \ xEj 



p / r>(p~<)( 



Qf } {E, z) = / dE p N p (E p , z) u pl {E p , z) 



is the distribution of pions produced per pp collision, with * ' E 
x = E/E p , = a' + 0.25, a' = 3.67 + 0.83£ + 0.075L 2 , x ^A^p) ~ °- 2 ^ P ) 

r = 2.6/v^ 7 , and a = 0.98/V^ 7 (see Kelner et al. 2006). = 5 N v$E, z) uJ pl {hE v , z) Af n (5E n ). (35) 
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Here the P7 collision frequency is for all charged pions is N^(E w ,z) — N Tr +(E 7r ,z) + 

^ 2e N w - (E n , z), it follows that 

c P7 (£ p ,z) = ±- 1°° de^^l [ ^ de'a${e>) e', (36) (m „, 

and the mean number of positive or negative pions is 

dn 

AfnP ,Pl + 2p2 , (37) >< N^E^z) dE ^ {E„E„). (44) 

with pi and p 2 = 1 — pi as the probabilities of single pion Similarly, 
and multi-pion production, respectively. These are related 

to the mean inelasticity function K n ^ — t~}oj~} by f E , 



71 7T,dec 



E u 



P1 = K " K *\ (38) dn 



x N n (E„z) ^ ^ (E^Ek). (45) 



where = 0.2 and if 2 = 0.6. 

For illustration, we show the obtained pion and muon 

4.2. Steady-state distribution of charged pions distributions at z = z Q in Fig. for the cases of production 

caused by pp and pj interactions. In these plots, we also 
The steady pion distribution obeys the transport equation show the particle distributions that correspond to no energy 



T6|) with the replacement i es J — > t^ 1 (E, z). The solution is losses. The solution corresponding to no energy losses will 

simply have the form 

N ' iE - z) = \bmC ' dE ' Q ^ w -» (E ^ = F#fy< (46 » 

xexp { - T „( E , E <». (39) £^0 
with <m (^ z ) 

, f E dE"t~ 1 (E z) We also note that the muon distributions shown include 

Tn(E , E) = I — — - — . (40) the contributions of the muons with different hclicity states 

Je ' \ b ^ h >\ added up. 

Depending on whether WG US6 Qtt (E,z) or Q ( r\E,z) m 

this last expression, we obtain N^ pp) (E, z) or nI pi) {E, z). 4.4. Neutrino emission 

The total emissivity of neutrinos, 

4.3. Muon steady state distribution ^ / \ ^ m \ , ^ m ^ 

Q u (E,z) = Qi T -> 1/ (E,z) + Q (1 ->v(E,z), 

As discussed in Lipari et al. (2007), to take the muon energy . . ,. , ., ? ,. • , , 

, . , , . / j ,. is the sum of the contribution of direct pion decays plus 

loss into account, it is necessary to consider the production ^ of de s 

of left handed and right handed muons separately, which a muon ecavs - 
have different decay spectra: 



"HEM = J*), \ 2 e(x - r T ) (41) 



Q^(E,z)= / dE n t-\ ec {E^)N v {E n ,z) 
Je 



dE^ v E*x(l-r*) s '■ y 71 J y ' 6(1 



-^f^(E»;E w )= [X r ' J 6(x-r,), (42) 

fl-E^ S ff x(l - r^y with a; = E/E n , and 

with x = E^/Ek and rv = (to^/to^) 2 . 4 „B max 

The injection function of negative left handed and pos- Q^„(E, z) = ^ / — -t" ^(E^N^ (E^z) 

itive right handed muons is i=1 J£ -^m 



(48) 



,.'5 



(in 



5 „ 2 4 , / 2 1 8ar , , 

3x 2 + -x 3 + 3x 2 h,; 

3 3 V 3 3 

{-,+} 



(49) 



x I N„-(E„,z) ""^ ^ r '+\ and 



In this last expression, x = E/E^, M{i.2} = A 1 ! ' ^{3.4} 



dE, 



ii 



dn ^ + ^ _ \ k {l-2} = -^{3,4} = "I, (50) 

+N V+ (E 7V , z) — Mr (E^Ek) j . (43) acC ording to Lipari et al. (2007). 

' The neutrino intensity (in units of GeV _1 s _1 ), 

Because CP invariance implies that dn^^^^- /dE^ = r 

dn^+^^+fdEfj,, and since the above distribution obtained ^{E) — J^d r Q v {E,z), (51) 



8 



M. M. Reynoso and G. E. Romero: Magnetic effects on //-production in MQs 



a = 100, z = z 



a = 100, z = z 



> 

CD 

(3 





Fig. 7. Pion and muon distributions at the base of the jet in the left and right panels respectively, originated by pp 
interactions (black lines) and by prf interactions (green lines). Solid lines: distributions obtained considering cooling. 
Dashed lines: distributions obtained neglecting cooling. 



is shown in Fig. [8] for the different values of a and the jet 
half-opening angle £. 

The differential flux of neutrinos arriving at the Earth 
can be obtained as 



d ®» - 1 T (F\ 



(52) 



This quantity, weighted by the squared energy, is shown in 
Fig. [U for a source at a distance d = 2 kpc, different values 
of the jet opening angle, and different values of a. As a 
guide, we also include a typical upper limit as derived from 
AMANDA-II data, as well as the expected sensitivity for 
the next generation neutrino telescope (Halzen 2006, see 
also Aiello et al. 2007). 

5. Neutrino production through jet-wind 
interactions 

In high-mass microquasars, the donor star can present a 
strong wind with dense clumps of matter (Romero et al. 
2007). In this section, we apply a very simplistic model to 
estimate the possible high-energy neutrino emission pro- 
duced by the interaction of such clumps with jet matter. 

We consider that the matter composing the clump that 
is successfully interacting with the jet is uniformly dis- 
tributed within a slice of the jet with a thickness equal 
to the radius of the clump. This enables us to apply the 
same method we used in the previous sections, where the 
one-zone acceleration region was placed in the vicinity of 
the compact object. In this case, the acceleration zone will 
be located at high distances from the compact object, say 
around half the distance to the companion star. There, the 
density of jet matter is two or three orders of magnitude less 
than that of the clumps, and the magnetic field is expected 
to be much weaker than close to the compact star. 



We assume that the clump density is p c = 10~ 12 g cm~ 3 , 
and we consider two different radii: R c = 10 11 cm and R c = 
10 10 cm. The relativistic particles are accelerated in a region 
of length R c . We assume a lower acceleration efficiency, r\ = 
0.01, and we can calculate, as above, the different cooling 
rates. In doing this, we consider that the density of cold 
matter in the acceleration zone is 



n c {z) = - 



4 it 2 / 



3 r?m p 



(53) 



For illustration we show in Fig. [10] the obtained accelera- 
tion and cooling rates, in the case of a successful jet-clump 
interaction at z = 10 12 cm. 

We consider that clumps can penetrate the jet at dis- 
tances > 5 x 10 11 cm, where the density of cold protons in 
the jet begins to decrease below 10 9 cm~ 3 . The protons and 
electrons are then injected in a slice of thickness R c using 
expression ([B]), and normalizing it through Eq. ([7]) using 
<he\ = 0.1 and a = 100. The steady state distributions of 
protons and electrons are found using expression (|18[) tak- 
ing T CBC = R c /v h . Next, the injection of pions is found to 
be dominated by the contribution of pp interactions, given 
by expression (|33|) . 

According to Fig. \TU[ the decay process of the pions and 
muons dominate in almost all the relevant energy range, so 
in a first approximation we can neglect losses so that 



N w (E,z) 



(E,z) 



t*\E,z) 



Substituting this pion distribution in expressions 
(1451). we can obtain 



%(E,z) 



(E,z) 



tp(E, z) 



(54) 



and 



(55) 
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Fig. 8. Different contributions of the neutrino intensity produced at the base of the jet. The cases of a = 100 and a = 1 
are shown in the top and bottom panels, respectively. Black lines correspond to £ = 5° and green lines to £ = 1.5°. 
Left panels: contributions from direct pion decays (solid lines) and from muon decays (dashed lines). Right panels: 
contributions due to pp interactions (solid lines) and due to pry interactions (dashed lines). 



Following the steps discussed in Sect. 14.41 the neutrino 
intensity and differential flux can be obtained. The cor- 
responding differential flux of neutrinos weighted by the 
squared energy is shown in Fig. 111! There we show the re- 
sults corresponding to the half-opening angles £ = 1.5° and 
£ = 5° for R c = 10 10 cm and R c = 10 n cm. A duty cycle of 
10% was adopted, meaning that the clumps interact with 
the jet successfully, on average 10% of the time. The neu- 
trinos produced are to be observed with cubic kilometer 
detectors over a long period of time (several years). 

6. Discussion 

The presence of an equipartition magnetic field in the jets 
of microquasars implies a strong attenuation in the high 
energy spectra of pions and muons that could be produced 



by hadronic interactions. This effect is found to be relevant 
in the vicinity of the compact object, as can be seen in Fig. 
[7] The neutrino flux expected in this case is significantly 
reduced at energies above 1 TeV, which is the range to be 
probed by upcoming neutrino telescopes such as IceCube 
(see Fig. El). 

Due to the lifetime of pions being shorter than that 
of muons, neutrinos produced by the direct decay of pi- 
ons dominate over those originated by muon decays, since 
muons lose a significant fraction of their energy by syn- 
chrotron radiation before decaying (see Fig. [51 left panels). 

With o=l, i.e. for equal power carried in rclativistic 
protons and electrons, the neutrino contribution due to pp 
interactions is dominant at low energies, E < 10 TeV for 
£ = 1.5° and E < 3 TeV for £ = 5°. At higher energies, 
the P7 contribution becomes important (see Fig. [HJ bot- 
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E [GeV] E [GeV] 

Fig. 9. Differential neutrino fluxes weighted by the squared energy. The cases of £ = 5° and £ = 1.5° are shown in the 
left and right panels, respectively. Black lines correspond to a — 100 and green lines to a = 1. Solid (dashed) lines: losses 
of secondary pions and muons considered (neglected). 



10" 7 




E [GeV] 

Fig. 11. Differential neutrino fluxes weighted by the 
squared energy. The solid lines correspond to £ = 5° and 
the dashed ones to £ = 1.5°. The black lines correspond to 
i?, c = 10 11 cm and the green lines to R c = 10 11 cm. 

torn right panel). However, in a case with a = 100, which 
seems more realistic (Heinz 2006), pp interactions provide 
the most relevant mechanism for neutrino production. 

For windy microquasars, an additional neutrino contri- 
bution can arise from jet-wind interactions. In this case, 
the clumps composing the stellar wind could interact with 
the jet at large distances from the compact object (z > 



5 x 10 11 cm), where the jet particle density is much lower 
than that of the clumps. The magnetic field in those regions 
of the jet is expected to be much lower than at the jet base. 
This leads to a negligible synchrotron energy loss of sec- 
ondary pions and muons, and probably to the production 
of a neutrino flux whose detectability depends on several 
factors such as the density and size of the clumps, and the 
duty cycle corresponding to this type of interactions. 

The jet half-opening angle is another parameter that 
matters. We considered two cases: £ = 5° and £ = 1.5°. 
The first value is often assumed in the literature and the 
second one is another possibility of a more collimated out- 
flow (notice that for SS433, £ « 0.6°). For wide opening 
angles, the magnetic energy density is lower than for nar- 
row ones, and the magnetic field is also lower. This leads to 
a lower synchrotron loss rate and hence to a higher max- 
imum energy of the particles. However, for wide opening 
angles, the density of cold protons is lower, and the spec- 
trum of produced secondary particles is lower. These effects 
can be seen in Figs. 151 and fTTT where the weighted fluxes for 
£ = 5° are lower and more slowly decreasing with energy 
than those for £ = 1.5°. 

It can also be noted from these plots that, for R c = 
10 11 cm, £ = 1.5°, and a duty cycle of 10%, the expected 
integrated neutrino signal above 1 TeV from jet-clump in- 
teractions will make the neutrino output dominate the pro- 
duction at the base of the jet. In the case of low-mass micro- 
quasars; however, only the latter will be present (Romero 

6 Vila 2008). 

7. Summary 

We have studied the effects caused by the magnetic field on 
the secondary pions and muons that could be produced in 
microquasar jets. First, assuming an equipartition magnetic 
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Fig. 10. Proton and electron accelerating and cooling rates (top panels), and pion and muon cooling rates (bottom 
panels), for jet-clump interactions at z = 10 12 . Plots shown: accelerating rate (solid lines, top panels), decay plus escape 
rate (solid lines, bottom panels), adiabatic cooling rates (dotted lines), synchrotron cooling rates (dashed lines), pp 
cooling rate (short-dashed line, top left panel), irp cooling rate (short-dashed line, bottom left panel), pj cooling rate 
(dash-dotted line, top left panel), ir-f cooling rate (dash-dotted line, bottom left panel), IC cooling rate (dash-dotted line, 
right panels) 



field, we calculated the neutrino production at the base of 
the jet adopting the one-zone approximation. In this case, 
protons and electrons are shock-accelerated in a localized 
region near the compact object. A fraction q re \ ~ 0.1 of the 
kinetic power of the jet is transferred there to relativistic 
particles. The energetic protons cool mainly by adiabatic 
expansion, pp and synchrotron radiation. Hot electrons cool 
mainly by synchrotron radiation. Secondary pions can be 



produced through pp and pj interactions, where the targets 
are the cold protons and the synchrotron photons in the jet. 
The pions produced also lose energy by adiabatic expan- 
sion, and mainly through synchrotron radiation. Pions still 
decay giving muons and neutrinos, but with an attenuated 
spectrum at high energies, due to the effect of pion syn- 
chrotron losses. Muons also cool significantly before decay- 
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ing, giving a much lower neutrino contribution than what 
was expected from previous, simpler calculations. 

Finally, we discussed the case of neutrino production 
through jet-wind interactions in high-mass microquasars 
with clumpy stellar winds. The clumps can successfully in- 
teract with the jet at large distances from the compact ob- 
ject (z > 5 x 10 n cm), where the magnetic field is much 
weaker than at the jet base. We applied a very simple model 
adopting the one-zone approximation at different distances 
along the jet where the clumps can interact. The size of 
the acceleration region is taken as the radius of the clumps 
R c w 10 10 — 10 11 cm, and a duty cycle of 10% is assumed. 
The produced pions and muons in this case do not undergo 
significant cooling because the magnetic field is relatively 
low. Hence, the neutrino spectrum obtained is not mod- 
ified by synchrotron radiation of the secondaries at large 
distances from the compact object. 

As a conclusion, we find that the main contribution to 
the neutrino emission stems from pp interactions. Pions and 
muons produced in the vicinity of the compact object are 
strongly affected by synchrotron losses, and their spectra 
are attenuated at high energies. As a consequence of that, 
we find that the neutrino flux is much less than expected 
when these effects are not taken into account. An additional 
neutrino contribution, arising from jet-wind interactions in 
high-mass microquasars, is not affected by these magnetic 
effects. Still, the detection of a neutrino signal from mi- 
croquasars seems difficult, but not impossible with next- 
generation neutrino telescopes such as IceCube, depending 
on the specific parameters of the system. 

The simple models presented here serve to illustrate 
these effects, which are crucial for assessing the detcctabil- 
ity of a neutrino signal from this type of sources. A more re- 
alistic treatment including the convection of particles in the 
jets and a consistent description of the acceleration mech- 
anism will be presented elsewhere. 
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